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Abstract 

Continuing work initiated in an earlier publication [Sato and Asada, PAS J, 61, L29 
(2009)], we consider light curves influenced by the orbital inclination and eccentricity 
of a companion in orbit around a transiting extrasolar planet (in a planet-satellite 
system or a hypothetical true binary). We show that the semimajor axis, eccentricity 
and inclination angle of a 'moon' orbit around the host planet can be determined by 
transit method alone. For this purpose, we present a formulation for the parameter 
determinations in a small-eccentricity approximation as well as in the exact form. As 
a result, the semimajor axis is expressed in terms of observables such as brightness 
changes, transit durations and intervals in light curves. We discuss also a narrow 
region of parameters that produce a mutual transit by an extrasolar satellite. 

Key words: techniques: photometric — eclipses — occultations — planets and 
satellites: general — stars: planetary systems 

1. Introduction 

It is of general interest to discover a second Earth-Moon system. Detections of such 
an extrasolar planet with a satellite (or hypothetical binary planet systems that do not exist 
in the Solar System) and probing the nature of such objects will bring important information 
to planet (and satellite) formation theory (e.g., Williams et al. 1997, Jewitt and Sheppard 
2005, Canup and Ward 2006, Jewitt and Haghighipour 2007). If a giant planet with a (perhaps 
Earth-size) rocky satellite were located at a certain distance from their host star, the satellite 
may be habitable and show vegetation, though these issues are out of the scope of this paper. 

It is not clear whether the lAU definition for planets in the Solar System can be applied 
to extrasolar planets as it is. The lAU definition in 2006 is as follows: A planet is a celestial 
body that (a) is in orbit around the Sun, (b) has sufficient mass so that it assumes a hydrostatic 
equilibrium (nearly round) shape, and (c) has cleared the neighborhood around its orbit. 

Regarding (c), the Earth can be called a planet, mostly because the common center of 



mass (COM) of the Earth-moon system is below the surface of the Earth. On the other hand, 
the COM of the Pluto-Charon system is located above the surfaces of these objects. Therefore, 
it is interesting to determine the COM position of a planet-companion system. In order to 
determine it, we have to know the true (not apparent) distance between the two objects. For 
this reason, Sato and Asada (2009) considered extrasolar mutual transits, as a complementary 
method of measuring not only the radii of two transiting objects but also their separation 
(See Sato and Asada 2009 also on detection probabilities of extrasolar mutual transits and a 
possible limit by Kepler Mission). As a particular case, a short separation binary, which has 
a rapidly orbiting companion, gives us a unique opportunity to measure the true separation of 
the binary, whereas a long separation one gives us only the apparent separation. Their work is 
very limited, however, in the sense that they assume circular orbits and also coplanar orbits as 
/ = 90 degrees for both planet and moon. Clearly it is important to take account of the orbital 
inclination and eccentricity. The main purpose of this paper is to study effects of the orbital 
inclination and eccentricity of a companion on extrasolar mutual transits. 

Since the first detection of a transiting extrasolar planet (Charbonneau et al. 2000), 
photometric techniques have been successful (e.g., Deming et al. 2005 for probing atmosphere, 
Ohta et al. 2005, Winn et al. 2005, Gaudi & Winn 2007, Narita et al. 2007, 2008 for measuring 
spin-orbit alignment angle). In addition to COROT^, Kepler^ is monitoring about 10^ stars with 
expected 20 ppm (= 2 x 10~^) photometric differential sensitivity for stars of V=12. This will 
marginally enable the detection of a moon-size object. In fact, COROT detected a transiting 
super-Earth (Leger et al. 2009, Queloz et al. 2009). 

Sartoretti and Schneider (1999) first suggested a photometric detection of extrasolar 
satellites. Cabrera and Schneider (2007) developed a method based on the imaging of a planet- 
companion as an unresolved system (but resolved from its host star) by using planet-companion 
mutual transits and mutual shadows. As an alternative method, timing offsets for a single 
eclipse have been investigated for eclipsing binary stars as a perturbation of transiting planets 
around the center of mass in the presence of the third body (Deeg et al. 1998, 2000, Doyle et al. 
2000). It has been recently extended toward detecting 'exomoons' (Szabo et al. 2006, Simon et 
al. 2007, Kipping 2009a, 2009b). Sato and Asada (2009) investigated effects of mutual transits 
by an extrasolar planet with a companion on light curves. In particular, they studied how the 
effects depend on the companion's orbital velocity. Furthermore, extrasolar mutual transits 
were discussed as a complementary method of measuring the system's parameters such as a 
planet-companion's separation and thereby of identifying them as a true binary, planet-satellite 
system or others. 

Their method has analogies in classical ones for eclipsing binaries (e.g., Binnendijk 1960, 
Aitken 1964). A major difference is that occultation of one faint object by the other transiting 

1 http: / /www.csa.int /SPECIALS /COROT / 
^ http://kepler.nasa.gov/ 
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a parent star causes an apparent increase in light curves, whereas echpsing binaries make a 
decrease. What is more important is that, in both cases where one faint object transits the 
other and vice versa, changes are made in the hght curves due to mutual transits even if no 
light emissions come from the faint objects. In a single transit, on the other hand, thermal 
emissions from a transiting object at lower temperature make a difference in light curves during 
the secondary eclipse, when the object moves behind a parent star as observed for instance for 
HD 209458b (Deming et al. 2005). 

This paper is organized as follows. In section 2, we consider effects of the orbital in- 
clination and eccentricity of a companion on light curves. For simplicity, henceforth, such a 
companion orbiting around a host planet is called a 'moon' even if it is not a satellite but a 
component of a hypothetical binary planet. In section 3, we present a formulation for param- 
eter determinations. Some numerical examples are also presented. Section 4 is devoted to the 
conclusion. 

2. Effects of the Orbital Eccentricity of a Transiting 'Moon' on Light Curves 

2.1. Approximations and notation 

The time duration of a transit, say a few hours, is much longer than the orbital period 
of an extrasolar planet, say a few days or greater. In one transit, the effect of the motion of 
the moon is much larger than that of the planet orbiting around a star. During the transit, 
therefore, we employ a constant velocity approximation only for the orbital motion of a planet- 
moon system around their host star. For a short separation case, on the other hand, we take 
account of the eccentric orbit of the moon, because the orbital period of such a moon around 
the planet may be comparable to (or shorter than) the timescale of the transit. 

The co-planar assumption that the orbital plane of a moon around its primary object 
is the same as that of the planet in orbit around the host star seems reasonable because it 
seems that planets are born from fragmentations of a single proto-stellar disk and thus their 
spins and orbital angular momentum are nearly parallel to the spin axis of the disk. Irregular 
satellites such as Triton, however, have significant inclinations presumably through capture 
processes. This requires that we should include the effect of orbital inclinations. We assume 
only the moon's orbital inclination, because it makes substantial effects on mutual transit 
light curves. Inclinations of the planet's orbital plane have been well understood and already 
observed (Charbonneau et al. 2000). 

Here we list our assumptions for clarity. 

• The inclination angle of the COM of the planet and moon is fixed at 90 degrees. 

• We take account of the inclination angle of the moon's orbit. 

• We assume that the planet-moon COM has a constant velocity (during a transit by the 
planet in front of the host star). 
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• The longitude of ascending node of the moon equals zero. 

• The planet-moon COM orbit has zero eccentricity. 

• We assume no limb darkening effects. 

For an eccentric orbit, we generally have the Kepler equation as 
T 

t = to + —{u-esmu), (1) 
ivr 

where to, T, e and u denote the time of periastron passage, orbital period, eccentricity and 
eccentric anomaly, respectively (e.g., Danby 1988, Roy 1988, Murray and Dermott 1999). In 
the following, we use the true anomaly / instead of the eccentric anomaly u (e.g., Danby 1988, 
Roy 1988, Murray and Dermott 1999). They are related by 




1 + e u 

tan-. (2) 

1-e 2 ^ ^ 

The distance between the orbiting body and a focus of the ellipse is written as 

a(l-e^ 

^"1 + ecos/- ^"^^ 
We use these equations for describing a moon orbiting around a planet. We denote the 
mean motion of the moon in orbit around the primary as rim = 211 /Tm, where the subscript m 
means the moon's quantity. The subscript p denotes the planet's quantity. 

For investigating transits, we need the transverse position x and velocity v of each object. 
We denote those of the COM for planet-moon systems as xcm and vcm, respectively, where the 
origin of x is chosen as the center of the star. We assume vcm as constant during the transit. 
The position and velocity of each planet with mass Mp and Mm in the planet-moon system as 
Xi and Vi {i = p,m), respectively. The direction of the observer's line of sight is specified by the 
argument of pericenter as an angle denoted by Um (See also Fig. 1). We express the transverse 
position as 

Xcm = Vcm it-tcu), (4) 



X. 



p 



Xcm + Op [(cos M - e„) cosw^ ^ 1 - e^siuMsinw^], (5) 



Xm = Xcm - am[{cosu - em) cosum. + V 1 - e^sinnsinwm], (6) 

where the semimajor axis of the orbit of each object around their COM is denoted by Oj, and 
tcM means the time when the binary's common center of mass passes in front of the center of 
the host star. In terms of the true anomaly, they are rewritten as 

, (l-e^)cos(/ + a;^) 
Xp — xcM-T-ap — , {I) 

{1 - el) COsjf + Um) /^x 

Xm = Xcm - am — z • (8) 

1 + emCos/ 

(See Table 1 for a list of parameters and their definition). 

The azimuthal velocity of the secondary object around the primary is 
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Vf = r— 
^ dt 



a n ^ + ^"""^^ (9) 
1 — e? 



m 



where / denotes the true anomaly (Murray and Dermott 2000). Here, apm and denote the 
semimajor axis and eccentricity of the moon's orbit with respect to the host planet. We assume 
that both the eccentricity of COM orbit vanishes and the inclination angle of COM equals 90 
degrees. Hence we can avoid a careful treatment of "sky-projected transverse position" . 

2.2. Transits in light curves 

We denote the intrinsic stellar luminosity as L. The apparent luminosity L' due to 

mutual transits is expressed as 

, , S-AS ,^ , 

L=Lx— (10) 

where 5" = vri?^, Sy = nRp, Sm = T^Rm^ = Sp + Sm — Spm- Here, Rs, Rp and Rm denote 
the radii of the host star, planet and moon, and Spm denotes the area of the apparent overlap 
between them, which is seen from the observer. Without loss of generality, we assume that the 
primary is larger than the secondary as Rp > Rm- 

2.3. Effects on light curves 

We investigate light curves by mutual transits due to planet-moon systems. The orbital 
velocity is of the order of apmUm- Therefore, we have two cases; vqm < o-pmnm and vcm > o.pm'n.m- 
The dimensionless ratio of the moon's orbital velocity to the planet's one is defined as 

W = "^""^"^ . (11) 
Vcm 

If Vcm < ttpmnm, we call it a fast case. If vcm > CLpmnm, we call it a slow one. The Earth- 
Moon (W = 0.03), Jupiter-Ganymede (W = 0.8) and Jupiter-Io (W = 1.3) systems represent 
slow, marginal and fast cases, respectively. Figure 2 shows a schematic light curve by mutual 
transits. 

Fig. 3 shows a slow case in circular motion, where we assume Rg'. Rp: Rm = 20 : 2 : 1. We 
assume also the same mass density for the two transiting objects and hence obtain Op : = 1 : 8. 

Eccentric orbit cases {W = 6 and Cm = 0.3) are shown by Figs. 4, 5 and 6 {um = vr/2, 
and 7r/4, respectively). Some parameters are chosen so that effects in the figures can be 
distinguished by eye, though such an event is unlikely to be detected by current observations as 
discussed later. For generating the ingress and egress of the various parts of the lightcurve, we 
do not use a linear interpolation but compute numerically the apparent overlap area between 
the objects. Here we assume the same configuration except for the observer's line of sight. For 
simplicity, we take to = tcM = in these figures. 

These figures show also the transverse positions of transiting objects with time, which 
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would help us to understand the chronological changes in the light curves. In particular, it can 
be understood that such characteristic patterns appear only when two objects are in front of 
the star and one of them transits (or occults) the other. 



3. Formulation for Parameter Determinations by Transit Method 

3.1. Parameter determinations from transit observations alone 

In all the above cases, the amount of decrease in light curves or the magnitude of 
fluctuations gives the ratios among the radii of the star and two faint objects {Rs,Rp,Rm)- The 
decrease ratios in the apparent brightness due to transits by the planet and moon are written 
as 



The stellar radius Rs (and mass Ms) are known for instance by its spectral type. Hence, the 
radii are expressed in terms of observables Rg, Ap and as 



We define the ratio between the brightness changes by the two objects as 



Circular Orbit and Orbital Inclination: 

First, we discuss a circular orbit in order to simply explain our idea. For more rigorous treatment 
of eccentric orbits, please see below, where we will finally give expressions for determining the 
separation apm- 

Behaviors of apparent light curves depend on W. Therefore, apmnm (as its ratio to 
vcm) can be obtained (Sato and Asada 2009). Seager and Mallen-Ornelas (2003) presents an 
analytic solution of parameter determinations for a single transit in the circular orbit case. 
Their solution can be used for our case of mutual transits by a planet and moon in front of a 
host star. In our case, their equations are rewritten as follows. The duration of the 'fiat part' 
of the transit (t^m) is described by 




(12) 



(13) 




(14) 
(15) 




(17) 



whereas the total transit duration (trm) is done by 




(18) 
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where Im denotes the orbital inchnation angle of the moon. By combining these equations, the 
impact parameter (bpm) can be derived as 



7 ^pm J 

'^pm — J-) COS 1 ^ 
Up 



A)2 - [sin^(t^^7r/T„)/sin^(iTm7r/T^)](l + VA)^ ^^^^ 



^ 1 - [sin^itFrnT^/Tm)/ Sin'^ (tTm'^/Tm)] 

The ratio apm/ Rp can be derived directly from Eq. (18) as 



(20) 



;i + v/A)2 - [1 - sin^(tr„^7r/T„ 

Rp \ sin^(tTmVr/Tm) 

Therefore, one can obtain apm as Rs x {Rp/ Rs) x (apm/Rp)- 

With ttpm. and rim in hand, one can thus estimate the total mass of the binary by 
GMtot = n'^ttp^ from Kepler's third law, where G denotes the gravitational constant. The 
orbital velocity apmnm gives the mutual force between the binary. 

If we assume also that the mass density is common for two objects constituting the 
binary (this may be reasonable especially for similar size objects as Rp ~ Rm), each mass is 
determined as Mp = Rp{Rp + Rl^)~^Mtot and Mm = Rm{Rl + Rm)~^ Mtot, respectively. Therefore, 
the orbital radius of each body around the COM is obtained as Op = R^{Rp + R%^)~^apm and 
ttm = Rp{Rl + R^m)~^(^pm, respectively. At this point, importantly, the two objects can be 
identified as a true binary (a^ > Rp) or planet-satellite system (a^ < Rp). However, a gaseous 
giant planet with a rocky satellite would exhibit largely different densities and this may be one 
of the most likely scenarios. 

In a slow spin case, on the other hand, the apparent separation a± (normal to our line 
of sight) is determined as a_i = Ti2Vcm from measuring the time lag T12 between the first and 
second transits because vcm is known above (Sato and Asada 2009). 
Eccentric Orbit and Edge-on Case: 

Henceforth, we take account of the orbital eccentricity of a moon for an edge-on case. In 
this case, intervals between neighboring "hills" are not constant because of the eccentricity. 
However, a time duration between three successive "hills" is nothing but the orbital period of 
the moon. Therefore, one can measure the period T^. We obtain Um as 

nm = (21) 

A key idea for determining the eccentricity is as follows. As for timescales, we have 
two observable ratios as T2/T1 and T12/T21. The former is the ratio between the widths of 
neighboring hills, whereas the latter is that between the transit intervals. On the other hand, 
we have two additional parameters and Um to be determined. Importantly, the number 
of measurable ratios is the same as that of the parameters that we wish to determine. In 
principle, therefore, the above two ratios may allow us to determine the two parameters Cm and 
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cjm, separately. This will be discussed in detail below. 

To be more precise, the full width of a "hill" at top and bottom are expressed as (See 
also Figure 2) 

T,op = ^^^V^' (22) 

_ 2{Rp + Rm) , . 

-I- bottom — -r-r ■ \^'-') 

Only for symmetric binaries {Rp = Rm), we have T^p = and thus true spikes. Otherwise, 
truncated spikes (or "hills") appear. 

For the primary transit, where the moon moves in front of the planet, we have / = 
TT/2 — Um- From Eqs. (9), (22) and (23), therefore, we obtain 



J itop — -— : , l^4J 

(^pm^m J- + Cm Sm UJm 



_ 2{Rp + Rrn) yi-e^ 

J-lbottom — T~ : • \^'^) 

^pva^m ^ ~r Cm Sm OJ^ 

For the circular orbit = 0, the second factors in the R.H.S. of these expressions become the 
unity and the first factors recover the case for = (See Sato and Asada 2009). 

For the secondary transit, where the moon moves behind the planet, we have / = 37r/2 — 
ojm- From Eqs. (9), (22) and (23), therefore, we obtain 



2{Rp — R 



1-e: 



T,top= ^' ^' -^ r^, (26) 



^pm'^m ^ Sill UJrn 



_ 2{Rp + Rrn) yl-g^ 
■J-2bottom — : • y^i) 

(Xprri^^m -L Sm LOfn 



We immediately obtain from Eqs. (24), (25), (26) and (27), 

rp _ ^2top 



(28) 



Itop 

T 

2bottom 

T 

Ibottom. 

Equation (28) is rewritten as 
Tr-1 

emSmUm = TfT—T 
J. J. ~\~ -L 

^Tr, (29) 

where the R.H.S. can be determined by observations alone. Once either or Um is known, 
Eq. (29) determines the other. 

Next, we consider the time interval between the primary and secondary transits. In 
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order to compute such an interval, one can use the Kepler's second law (the constant areal 
velocity). After lengthy but straightforward calculations, the area swept from the primary 
transit (/ = 7r/2 — Um) till the secondary (/ = 3n/2 — ujm) becomes 



S 



12 



TT 1 

— h arcsin if H — sin(2 arcsin ) 
2 2 ^ ' 



where 6^ denotes the semiminor axis of the moon's elliptic orbit and H is defined as 
if = 



COt^ UJr. 



l-e^ + cot^o;^' 
By using Eq. (29) in Eq. (31) for eliminating coto;™,, we obtain 



H 



R 



\ 1 



(30) 



(31) 



(32) 



It is convenient to use T12/T instead of T12/T21, because T12/T gives a simpler expression 
than T12/T21, which is used in the above explanation of the key idea. The total area is S" = 



irambm- Therefore, we find 



T12 
T 



— I — arcsmii H — H \ 
2 TT n ] 



(33) 



This includes without Um- Hence by using this relation, one can determine separately the 
eccentricity by measuring the time intervals. We should note that Eq. (33) is valid even for a 
general case with a certain inclination angle. 

Here, we consider a small eccentricity approximation, which may be useful for a quicker 
estimation of the parameters. For small Cm, Eq. (31) is expanded as 

H = err,cosujrn + 0{el^). (34) 

Substitution of this into Eq. (33) gives us 

Tl2 

T 



1 2^771 

- H cosoj^ 

2 TT 



0{el 



(35) 



The correction to the circular case (Ti2/Tm = 1/2) is 2emCosc<;.m/7r ~ O.Tcm coswm. Even for 
Um ~ for instance, it leads to seven percents for = 0.1. 
We define the difference between T12 and T21 as 



6T = T12 -T21. 

Replacement Um by Um + tt changes Eq. (35) into 
721 _ 1 _ 2e, 
T 



■cosUm + 0{el 



2 TT 

By using Eqs. (35) and (37), we obtain 



(36) 



(37) 
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ttST ^/ c{, 

emCosuJm = ^Tf^ + Oie-"). (38) 



' m 



Let us consider every cases separately. 

(I) If and only if the L.H.S. of Eqs. (29) and (38) vanish, Cm does. This means a circular 
motion and thus the observer's direction ujm becomes meaningless. 

(II) For a case when the L.H.S. of Eq. (29) vanishes but that of Eq. (38) does not, we 
find Cm 7^ and sinwm = 0, namely 

= (modvr). (39) 

Eq. (38) immediately gives 

Tx5T , , , 

^^ = T^- (40) 

Hence, the orbital eccentricity is determined. 

(HI) If the L.H.S. in Eq. (38) vanishes but that in (29) does not, we find 7^ and 

cosWm = 0, namely 
vr 

= - (modTr). (41) 

Eq. (29) immediately gives 

em = Tr. (42) 

Hence, the orbital eccentricity is measured. 

(VI) A general case in which the L.H.S. of neither Eqs. (29) nor (38) vanish: By dividing 

Eq. (38) by Eq. (29), we obtain 

Ti 5T 1 , , 

cota;^ = - — — , (43) 

where the R.H.S. can be determined by observations alone and hence this equation gives us the 
observer's direction Um- By substituting the determined Um into Eq. (29), one can find the 
value of the eccentricity. 

Up to this point, Cm and ujm both are determined. Eqs. (24) and (25) are rewritten as 



^ _ Tm{Rp — Rm} 

Tm{Rp — Rm) (Titop + T2top) V 1 ~ 6: 



o2 



(yRp Rm) ( TjYi Tjj 



27r y Titop T2top ^ 



+ Jl-el, (44) 



Tm{Rp + Rm) \/^ 

ttr. - 



"pm 



T^Tibottom 1 + Cm siu Um 



Tm{Rp + Rm) {Tlbottom + T2bottom)y 1 ~ 6: 
2 TtTi bottom T2bottom 
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p2 



(^Rp ~\~ Rm) ( T^m 



^71 \ i Ibottom 2bottom / 

respectively, where we used Eq. (29) If and only if Rp = R^, we obtain T^op = T2top = and 
Eq. (44) thus becomes undetermined, whereas Eq. (45) is still well-defined. 

When one wishes to consider the secondary transit instead of the first, one can use 

Tm{Rp — Rm) \/l~^m , . 

apm = 7^ : , (46) 

'^'L2top J-— CmSmWm 



o2 
■'m 



(47) 



TTi 2bottom J- — Cm Sm Um 

They are obtained by replacing Um with Wm + Tr in Eqs. (44) and (45). By noting Titopil + Tn) = 
T2bottom{X-TR) and T2top(l +T/j) = T2feottom,(l -7r), one can show that Eqs. (46) and (47) agree 
with Eqs. (44) and (45), respectively. By using one of these expressions, we can thus measure 
the semimajor axis of the eccentric orbit. In terms of the decrease in apparent brightness, Eqs. 
(44) and (45) are written as 



a. 



' \ ' (^ + ^]^l-el, (48) 

-Tts ^TT \I Hop 1 2top , 



VA^) / Try 



+ h/l-^L (49) 



Rs 27r \Tibottom T2bottom- 

where we used Eqs. (14) and (15). 

Determination of the semimajor axis is sensitive to measurement errors in the widths 
of the hills. This statement can be proven by using Eq. (49). For simplicity, we assume 
Ap ~ A™ ~ A and Tuottom ~ T2bottom ~ Tbottom, so that Eq. (49) can be reduced to 
2„ T„ 



apm ~ -RsVAjt^^Ji - el. (50) 

^ bottom 

The logarithmic derivative of this becomes 

dctpm dRg ^ c/A ^ dT^n dTbottom ^mdCm (51) 
Opm Rs 2 A Tm Tbottom (1 ^m) 

We focus on dTm/Tm and dTbottom. /Tbottom, because we can expect much more accurate mea- 
surements for Rs and A, and the last term involving may be relatively small (e^ < 1 and 
dcm < 1)- Time resolution in observations seems dTm ~ dTbottom, while ^ Tbottom- Therefore, 
dTm/Tm <^ dTbottom /Tbottom, which means that accurate measurements of Tbottom are crucial for 
the determination of apm- 

Figure 7 shows a flow chart of the parameter determinations that are discussed above. 
The above formulation for parameter determinations actually recovers the correct values in 
Figs. 4-6. In the numerical examples, the original parameters are retrieved within twenty 
per cents. 

Eccentric Orbit and Orbital Inclination: 
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Figure 8 shows a difference between liglit curves for tlie edge-on (/„ = 90 deg.) and an inclination 
case {Im = 88 deg.). Let 2;-axis denote tlie axis normal to the x-axis on the celestial sphere. We 
define the distance of a 'moon' from the 2-axis at the initial time of the mutual transit as 



Sb 



\ 



{Rp + R„ 



COS^ /„ 



(52) 



\ l + e^cos/^ 

where the subscript b means that the quantity is related with Tibouom and T2bottom as shown 
below (See also Fig. 9). Similarly, when the 'fiat part' of the spike in light curves starts (or 
ends), we define the distance of a moon from the 2;-axis at this epoch as 



St 



\ 



(^Rp Rn 



COS^ /„ 



(53) 



1 + e,„cos/ 

where the subscript t means that the quantity is related with Tuop and T2top as shown below. 
Therefore, we obtain the duration Ttop as 
2st 



T — i 

top -r - 1 



f 



where Vf is given by Eq. (9). For the primary transit (/ = 7r/2 — Um), we thus obtain 



2s 



ti 



Itop 



whereas for the secondary (/ = 37r/2 — w^), we have 



T, 



2s 



t2 



2top 



'2pm'^m(l Cm.sinCc'^ 

Here, we define Sa and St2 as 



Stl 



St2 



\ 



{Rp — R„ 



1 + ernsinur, 



COS^ /„ 



{Rp — R„ 



'2pm(l ^m) 



COS^ /„ 



\1 — Cm sincu^^ 

In the similar manner, we obtain the width of the spikes at the bottom as 

T 

J- hntf.nm. , ~ • 



f 



For the primary transit (/ = 7r/2 — Um), we thus obtain 



2s 



61 



Ibottom 



^pmX^m ( siu CU-^ 



whereas for the secondary (/ = 37r/2 — uom), we have 



2bottom 



Opr?i'^r?i(l Cm sin J] 



(54) 

(55) 

(56) 

(57) 
(58) 

(59) 
(60) 
(61) 



where we define Sbi and Sb2 as 
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Sbl 



\ 




{R, + Rmr- """^^^ ."-^ ) cos2 J^, (62) 



Sb2 



\ 



{R, + RmY- (^EI!^LJ^\ ' cos2 1^. (63) 

Because of the orbital inclination, we have to consider T^p and Tbouom, separately. We 
define the ratios as 

Trtop^^, (64) 

and 

T 

Ti ''-2bottom 
rbottom = • Vt)>3j 

Ibottom 

Substitutions of Eqs. (55), (56), (60), (61) into these ratios lead to 
_ St2 l + emsinum 

-Lrtop — — :; : ; l,t)t)J 

Stl 1 - CmSmUm 

_ S62 l + e,„sina;,„ 

J- rbottom — ' 1 • 1,0' j 

Sbl l-emsmum 

For the edge-on case (Jm = 90 deg.), we obtain St2/sti = Sb2/sbi = 1. Then, we have 
Trtop = Trbottom- For a general case (/^ 7^ 90 deg.), on the other hand, we find T^op 7^ ^bottom- 
We thus expect that a ratio between them will give us the information about the orbital 
inclination. The ratio is 



Trtop _ St2 Sbl 



(68) 



Trbottom Stl Sb2 

The L.H.S. can be measured by observations. 

There are four unknown quantities apm, Cm, hn and Um- We have four equations of (66), 
(67), (68) and the last one that can be chosen out of (55), (56), (60) and (61). Therefore, one 
can determine the quantities a^m, e^, Im and ujm by using these equations for observations. 
For practical observations, data fittings at the slope of light curves are used, instead of tran- 
sit durations, for determinations of the orbital inclination angle (Charbonneau et al. 2000). 
Nevertheless, an analytic solution is necessarily worthwhile to understand the properties of a 
given physical system, even when numerical fits are in practice the best way to determine the 
system parameters. 

A partial transit occurs if the apparent impact parameter of the moon is in {Rp — Rm.,Rp + 
Rm)- For the primary transit (/ = 7r/2 — Um), it occurs if the orbital inclination angle satisfies 

Rp — Rm 1 + Cm sinWm -Rp + Rm 1 + Cm siu /„„,, 

— YZr^ ^ cosim < T^-^ • (69) 

CLpm J- (^m dpm J- 

For the secondary one (/ = 3-7r/2 — Um), the condition of a partial transit becomes 

Rp Rm 1 (^m simUm y Rp ~^ Rm 1 ^m siu.UJm /i-rn\ 
— TZ~^ ^ COsIrn < — TZ~^ ■ ^^^^ 

CLpm J- (^m dpm J- (^m 
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Such a partial transit by a moon orbiting a host planet produces a 'U'-shaped spike in light 
curves (See Fig. 10). 

3.2. Timescales and brightness changes 

We have presented a formalism for parameter determinations. Before closing this section, 
let us make brief comments on typical timescales and amplitudes in the brightness changes. 
The timescale of a brightness change due to a giant planet is about 



Therefore, detections of such fluctuations due to mutual transits of extrasolar planet-moon 
systems require frequent observations, say every hour. Furthermore, higher frequency (e.g., 
every ten minutes) is necessary for parameter estimations of the system. 

Let us mention a connection of the present result with space telescopes in operation. 
Decrease in apparent luminosity due to the secondary planet is 0(i?^/i?^). Besides the time 
resolution (or observation frequency) and mission lifetimes, detection limits by COROT with 
the achieved accuracy of photometric measurements (700 ppm in one hour) could put Rm/Rs ~ 
2 X 10^^. The nominal integration time is 32 sec. but co-added over 8.5 min. except for 
1000 selected targets for which the nominal sampling is preserved. By the Kepler mission 
with expected 20 ppm differential sensitivity for solar-like stars with my = 12, the lower limit 
will be reduced to Rm/Rs ~ 4 x 10^^. An analogy of the Earth- Moon {Rm/Rs ~ 2.5 x 10~^, 
W ~ 0.03) and Jupiter-Ganymede {Rm/Rs ~ 4 x 10"^, W ~ 0.8) will be marginally detectable. 
Observations both with high frequency (at least during the time of transits) and with good 
photometric sensitivity are desired for future detections of mutual transits. 

3.3. Dynamical limit and constraints on W 

For Roche limit, we have 



where /3 denotes a numerical coefficient < (3 < 1 and Rh is Hill radius (See Domingos et 
al. 2006 for more detailed stability arguments by numerical computations). For simplicity, we 
assume Mg ^ Mp ^ Mm- Then, we have ttpm ~ dm and Rh is approximated as 



where Mg denotes a host star mass and dp denotes the orbital radius of a planet orbiting the 
star. Kepler's third law gives apm^m and vcm as 




(71) 




(72) 




(73) 




(74) 



(75) 
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Combining these relations, therefore, we find 

31/6^-1/2 (76) 



Ms 

If one assumes Mp ~ Mj (Jupiter mass), we obtain W > 10^^. This lower bound is less severe. 
On the other hand, there is a stringent constraint that the moon's closest approach r^m cannot 
be within the planetary radius. We thus have 

am > • (77) 

This leads to 



mm 



If one assumes the jovian mass and radius (Rj), this is rewritten as 



ylAUj \MJ \Mpy .,™„, 
For W = 6 and = 0.3, we obtain dp > 25 AU. Namely, a planet with a long orbital period 
Tp > 125 years is required. Kepler mission for several years is unlikely to see a transit by such 
a long period planet. 

Next, we consider a constraint that a mutual transit can occur. The transit duration for 
a planet in circular orbit is 

D^^, (80) 
ndp 

where we assume the maximum duration by taking the vanishing impact parameter (See Seager 
and Mallen-Ornelas for a more accurate form). As this limiting case for the Roche limit, we 
obtain the fastest case of a 'moon' as 
/ M \ 

(si) (^1) 

Using the Kepler's third law, this leads to 

Tm-[jj Tp. (82) 
For our fast case, we require < D, which gives a bound on /3 as 

^<3V»(A) . (S3) 
We substitute this into /3 of Eq. (76) so that we can obtain 



w> 



/ VM, 



\dj) \rJ \Mj) \Ms 
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lyear/ \ Rs J \MjJ \ Ms 
where dj denotes the mean orbital radius of the Jupiter. 

On the other hand, the dynamical arguments put an upper bound by Eq. (79). This is 
rewritten as 

Therefore, there exists a narrow band as shown by Fig. 11. Outside this range, the proposed 
method cannot work. 

4. Conclusion 

We have shown that light curves by mutual transits of an extrasolar planet with a 
'moon' depend on the moon's orbital eccentricity, especially for small separation (fast) cases, in 
which occultation of one faint object by the other transiting a parent star causes an apparent 
increase in light curves and such characteristic fluctuations with the same height repeatedly 
appear. We have also presented a formulation for determining the parameters such as the 
orbital eccentricity, inclination, semimajor axis and the direction of the observer's line of sight. 
This will be useful for probing the nature of the transiting planet-moon system. 

When actual light curves are analyzed, we should incorporate (1) photometric corrections 
such as limb darkenings, and (2) perturbations as three (or more)-body gravitating interactions 
(e.g., Danby 1988, Murray and Dermott 2000). 

We would like to thank the referee for invaluable comments on the manuscript. We 
would like to thank S. Ida, S. Inutsuka. E. Kokubo and Y. Suto for stimulating conversations 
and encouragements. This work was supported in part (H.A) by a Japanese Grant-in-Aid for 
Scientific Research from the Ministry of Education, No. 21540252. 
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Table 1. List of quantities characterizing a system in this paper. 



Symbol 


Definition 


T 


Orbital period of an extrasolar 'moon' around a planet 




Mean motion of a moon (= 21: /Tm) 


ClfYi 


Semimajor axis of a moon's orbit w.r.t. a planet-moon's center of mass 


dp 


Semimajor axis of a planet's orbit w.r.t. a planet-moon's center of mass 


CLpui 


Semimajor axis of a moon's orbit around a planet (= Op + 


a± 


Apparent separation of a planet-moon system 


Rs 


Radius of a host star 


Rp 


Radius of a planet 


Rm 


Radius of a moon 


Mp 


Mass of a planet 


Mm 


Mass of a moon 


Mtot 




XCM 


Transverse position of a planet-moon's center of mass 


VCM 


Transverse velocity of a planet-moon's center of mass 


Xp 


Transverse position of a planet 




Transverse position of a moon 


e 


Orbital eccentricity of the moon 


to 


Time of periastron passage of the moon 




Argument of pericenter of the moon 


tcM 


Time when the planet-moon's center of mass passes across the star's center 


Ap 


Decrease rate in apparent brightness due to the planet transit 


Am 


Decrease rate in apparent brightness due to the moon transit 


Tit op 


Time duration: width of a hill's top at the primary transit in light curves 


T\hottorn 


Time duration: width of a hill's bottom at the primary transit in light curves 


T2top 


Time duration: width of a hill's top at the secondary transit in light curves 


T2botto'm 


Time duration: width of a hill's bottom at the secondary transit in light curves 


Tl2 


Time lag from the first transit to the secondary 


T21 


Time lag from the secondary transit to the first 



5T T12 - T21 
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Fig. 1. Direction of the line of sight. It is denoted as Um, which is the argument of pericenter. 
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Tl bottom T2bottom Time 



Fig. 2. Schematic figure of a liglit curve due to a mutually transiting planet and moon in front of their 
host star. 
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W=0,W=1 
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98 - 
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-2 -1.5 -1 -0.5 0.5 1 1.5 2 

t 

Fig. 3. Light curves: Solid red one denotes tire zero limit of a 'moon' orbital motion as a reference 
(W = apmnm/vcM = 0). Dashed green one is a marginal spin case (large separation) for W = 1 (Sato 
and Asada, 2009). The vertical axis denotes the apparent luminosity (in pcrcents). The horizontal one 
is time in units of the half crossing time of the star by the COM of the binary, defined as Rs/vcm- If 
one takes Ms = Mq, Rg — Rq and vcm = 30 km/s (namely, 1 AU distance from the host star), t = 1 
corresponds to w 6 hours. For simplicity, we assume the binary with a common mass density, a radius 
ratio as Rs ■ Rp Rm = 20 : 2 : 1, and ap„i/Rs = 0.9. 
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W=6, e_m=0.3, omega=90 deg. 

100 I ' i ' ' ' ' 



99.5 - 




98.5 - 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

t 



W=6, e_m=0.3, omega=90 deg. 




t 



Fig. 4. Top panel: a light curve for a fast case (small separation) as = 6 and e,„ = 0.3. The radius and 
separation are the same as those in Fig. 3. The observer's direction is a;,„ = 7r/2. Brightness fluctuations 
appear with Titop/Tm = 0.013, Tuottom/Tm = 0.022, T2top/Tm = 0.039, T2bottom/Tm = 0.067, Ti2/T„ = 0.51 
and T2i/Tm = 0.49 (normalized by T,„, the orbital period of the moon: Tm =Ti2 + T'2i)- We obtain Tj. ^ 1.7 
from Eq. (29) and thus e„iSina;„i ^ 0.3, whereas Eq. (38) approximately gives us em COSWm ~ 0. Therefore, 
we recover well the parameters as ~ 0.3 and w„i ~ 7r/2, even in the linear approximation in e„i. Finally, 
Eq. (49) tells ttpm/Rs ^ 0.9. Bottom panel: motion of each body in the direction of x normalized by Ra 
(sohd red for the primary and dotted green for the secondary). When one faint object transits or occults 
the other in front of the host star, mutual transits occur and a 'hill' appears in the light curve. 
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W=6, e_m=0.3, omega=0 deg. 



100 



99.5 



98.5 




W=6, e_ni=0.3, oniega=0 deg. 




Fig. 5. Top panel: a light curve for all the parameters same as those in Fig. 4, except for the ob- 
server's direction as LOm ~ 0. Brightness fluctuations appear with Tuop/Tm = 0.018, Tibottom/Tm = 0.018, 
T2top/Tm = 0.050, T2botto,n/T„, ^ 0.050, Ti2lT^ = 0.69 and Tsi/T^ = 0.31. We obtain - 1.0 from Eq. 
(29) and thus emsinwm ~ 0, whereas Eq. (38) approximately gives us emCOSWm ^ 0.3. Therefore, we 
recover well the parameters as Cm ^ 0.3 and lo„-i ^ 0, even in the linear approximation in e,„. Finally, Eq. 
(49) tells apm/Rs ^ 0.9. Bottom panel: motion of each body (solid red for the primary and dotted green 
for the secondary). 
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98.5 



W=6, e_p=0.3, omega=45 deg 




-2 -1.5 -1 -0.5 



W=6, e_m=0.3, omega=45 deg. 




Fig. 6. Top panel: a light curve for all the parameters same as those in Fig. 4, except for the ob- 
server's direction as ujm = 7r/4. Brightness fluctuations appear with Tuop/Tm = 0.015, Tibottom/Tm ~ 0.022. 
T2top/Tm = 0.041, T2bottorn/T,n - 0.062, Tis/T^ = 0.62 and Tsi/Tm = 0.38. We obtain T, ^ 1.5 from Eq. 
(29) and thus em sinwm ^ 0.2, whereas Eq. (38) approximately gives us emcoswm ^ 0.2. Therefore, we 
recover well the parameters as e„i ^ 0.3 and Ltj„i ^ 7r/4, even in the linear approximation in e,„. Finally, 
Eq. (49) tells Upm/Rs ~ 0.9. Bottom panel: motion of each body (solid red for the primary and dotted 
green for the secondary). 
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Brightness changes 




+ R: 



Rp , R 



m 




Titop , Tztop 

+ (Tlbottom , Tzbottom) 

TiZjTzi 



+ R: 



Semimajor axis 



Fig. 7. Flow chart of parameter determinations. Starting from measurements of brightness changes, the 
semimajor axis ap„j can be finally determined for a fast case. 
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Fig. 8. Difference in light curves due to orbital inclinations of a 'moon'. All parameters but for the 
inclination angle are the same as those in Fig. 4. The solid (red) curve denotes /,„ = 90 degree case, while 
the dashed (green) one means a case of /,„ = 88 deg. Because of the orbital inclination, the duration of a 
mutual transit by the satellite is shortened. 
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Fig. 9. Definition of Sf, and St'. Tlie total transit duration Tbottom is given by Sb, wfiereas tlie 'flat part' 
of the hill Ttop is done by Sf. 
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Fig. 10. Light curve by a partial mutual transit of a satellite. All parameters but for the inclination 
angle are the same as those in Fig. 4. The solid (red) curve denotes /„, = 90 degree case, while the dashed 
(green) one means a ease of /,„ = 86 degree for a partial transit. Such a partial transit case produces a 
'U'-shaped hill in light curves. 
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Fig. 11. Possible bound on W and Tp. Here, wo assume Mg = A/©, Rg ~ Rq, Mp ~ Mj, rmin = Rj and 
e,n — 0. The sliaded (green) region denotes prohibited regions of the parameters that are constrained by 
Eqs. (84) and (85). 
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